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^ ^ We report the status of nucleon structure calculations on the (2+l)-flavor dynamical domain- 

qj wall fermions ensembles with pion masses as low as 180 and 250 MeV on a lattice with about 

f^h 4.6 fm spatial extent. A combination of the Iwasaki+dislocation- suppressing-determinant-ratio 

(I+DSDR) gauge action and DWF fermion action allows us to generate these ensembles at cutoff 
of about 1.4 GeV while keeping the residual mass small. Nucleon source Gaussian smearing has 
been optimized. Preliminary nucleon mass estimates are 0.98 and 1.05 GeV. 

00 
CO 



The XXVIII International Symposium on Lattice Field Theory, Lattice2010 
June 14-19, 2010 
Villasimius, Italy 

* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



2+lfDWF Nucleonfrom RBC/UKQCD 



Shigemi Ohta 



1. Introduction 

RIKEN-BNL-Columbia (RBC) and UKQCD collaborations have been investigating nucleon 
structure by calculating isovector form factors and some low moments of isovector structure func- 
tions [1, 2, 3, 4, 5, 6, 7]. 

In our 2+1 flavor dynamical domain- wall fermions (DWF) ensembles with the lattice cutoff 
of a- 1 = 1.73(2) GeV and spatial volume of (2.7fm) 3 [8], we found surprisingly strong volume 
dependence in isovector axial vector-current form factors [5, 6]: lattice spatial extent of five pion 
Compton wave lengths or larger seems necessary to reproduce these quantities. This provided the 
first concrete evidence of virtual pion cloud surrounding nucleon. While the vector-current form 
factors do not show such strong volume dependence, they fail to reproduce the experimental values 
as the pion mass in these calculations stayed much higher, at about 330 MeV, than the experimental 
value of about 140. It would be very interesting to do similar lattice-QCD calculations at lower 
pion mass and large spatial volume. 

On the other hand for low moments of isovector structure functions such as quark momentum 
or helicity fractions, our calculation at the lightest pion mass of about 330 MeV showed encour- 
aging trending toward the experimental values [7]: calculations at lighter quark mass values are 
necessary to confirm if this interesting trend is a real physics effect. 

Fortunately, the new 2+1 flavor dynamical DWF ensembles we started to generate last year 
at the lower cutoff of a~ l = 1.368(7) GeV [9] are ideal for such a study of nucleon structure at 
light quark mass and large spatial volume: the pion mass is about 180 and 250 MeV for the two 
ensembles and the spatial extent with 32 lattice units is physically about 4.6 fm. Here we report the 
current status of the nucleon structure calculations with these ensembles. 

2. Form factors 

The isovector form factors are conveniently summarized in neutron /3 -decay transition matrix 
elements, 

(p\V+(x)\n) =u p 

(p\A+(x)\n) = u p [y^FAiq^+iq^Fpiq 2 )] u n e iqx . (2.2) 

The vector-current form factors, the vector, Fy = F\, and tensor, Fj = Fi, are identical to the 
isovector components of the nucleon electromagnetic form factors, and are sometimes represented 

as electric, G^(g 2 ) = F\ =-F2, and magnetic, Gm = F\ + F2, form factors. From these we 

deduce such quantities as mean-squared charge radius and magnetic moment. These form factors 
have important implications in atomic physics as they determine nuclei interaction with photon 
and electron. They are under very active experimental studies: indeed a recent experiment [16] 
is in systematic disagreement with earlier ones about proton mean-squared charge radius. Our 
lattice-QCD numerical calculations of these quantities may help resolve this discrepancy. 

The vector form factor at zero momentum transfer is also related to the Fermi constant and 
Cabibbo mixing: gy = Fy(0) = GpermiCOsOcabibbo- The axial form factor at zero momentum trans- 
fer, the axial charge deviates from unity in units of gy because of QCD correction, g& = Fa (0) = 
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1.2694(28)gy [17]. Whether the lattice-QCD numerical calculation can reproduce this quantity is 
an interesting test of the method. The vector, gy, and axial, g&, charges determine the /3 -decay 
life time of neutron. The latter also dominates pion-nucleon interaction through the Goldberger- 
Treiman relation, mjygA x fngnNN- Thus these quantities determine primordial and stellar nucle- 
osyntheses and are important in explaining abundance of elements. 

To calculate these form factors in lattice QCD, we usen ratio of two- and three-point correlators 

. ^3pt ( f sink)0 . 

such as — - — -; — with 

C2pt(?sink) 

C 2p t(?sink) = £ ( {Np(t s - mk )N a (0)), (2.3) 
a.p V 1 J afi 



qpf (W) = £ r a p(N p (t sink )O(t)N a (0)), (2.4) 
a,j3 

where N denotes an appropriate nucleon operator, O the relevant current operator, and T an appro- 
priate spin or momentum-transfer projections. In this study we choose a standard nucleon operator, 
N = £abc{u T a CY5db)u c . For details of the spin or momentum- transfer projectors, T, we refer our 
earlier publications such as ref. [4]. The ratio gives a plateau in < t < t s - mk for a lattice bare value 
(O) for the relevant observable. We will discuss later how to optimize the sink position in time, 
4ink> to minimize systematics from excited-state contamination. 

Most recently we calculated all the four form factors using our 2+1 flavor dynamical DWF 
ensembles at lattice cutoff of a~ l = 1.73(2) GeV and spatial extent of 2.7 fm [5, 6]: at the lightest 
pion mass of about 330 MeV we found huge finite-volume corrections in the axialvector-current 
form factors. While the vector-current form factors do not suffer from such a finite-volume cor- 
rection, their values generally deviate from the experiments very likely because our degenerate up 
and down quark mass is still far from reality. Thus we like to calculate these quantities at more 
realistically light mass and at the same time in sufficiently large spatial volume. 

3. Moments of structure functions 



The structure functions are measured in lepton deep-inelastic scattering off nucleon. The 

a 2 

cross section is factorized in terms of leptonic and hadronic tensors, <=c ^,^ l^ v W uv , and because 

W ) 

we know the leptonic sufficiently well, = 2(k fl k' v +k' ll k v — -Q g^ v ), we deduce the hadronic 
tensor, W^ v ', from which we extract the structure functions: spin-unpolarized structure functions, 
F\ and Fi, are obtained from symmetrized hadronic tensor, 

w^H*,Q 2 ) = (-^ v + ^f) n(x,Q 2 ) + (> - (> - F ^f 1 , (3-D 

and spin-polarized ones, gi and g2, are from antisymmetrized tensor, 

W^\x,Q 2 ) = ieWq p ^( gl (x,«2 2 )+g 2 (x,«2 2 ))-^g 2 (^e 2 )) , (3.2) 
with v=q-P,S 2 = —M 2 , x = Q 2 /2v, and Q 2 = \q 2 \. 
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What are accessible for lattice-QCD calculations are moments of the structure functions, such 

as 

2 [ 1 dxx n - l F l (x,Q 2 )= £ c\ q) n {ll 2 /Q\g{mx n ) q {v) + 0{\/Q 2 ), (3.3) 

J ° q=u,d 

where c\ fl are perturbative Wilson coefficients in the operator product expansion, and {x n ) q (\L) are 
the moments we seek. The other unpolarized structure function, F% y is similarly expressed in terms 
of the same {x n ) q (li) moments but with different Wilson coeeficients C2, n - From polarized structure 
functions gi and g2 we obtain {x n )^ cl {pL) and d n {pL) moments. In addition to these moments of 
DIS structure functions, we calculate tensor charge, (l)s q = {P,S\ YiY5<3nvY\P,S) which may be 
measured by polarized Drell-Yan and RHIC Spin experiments. 

As the operator product expansion indicates, these moments need be renormalized before we 
can compare them with experiments. Good continuum-like flavor and chiral symmetries of domain- 
wall fermions makes this non-perturbative renormalization procedure feasible. We follow Rome- 
Southampton procedure [10, 11]. 

These moments can be calculated in lattice QCD as forward matrix elements of certain local 



operators, such as q 



74 £>4 -^£7*0* 

3 k 



q. Such matrix elements are obtained from lattice QCD 



as ratios of two- and three -point correlators just like in form factor calculations. We calculate those 
moments that do not require finite momentum transfer, namely (x) q , (1)a<?, {x)&q, and d\. For 
details we refer our earlier publication [4, 7] and references cited therein. 

We calculated these low moments using the same 1.7-GeV, 2.7-fm, (2+l)-flavor dynamical 
DWF ensembles [7]. In contrast to the axial vector-current form factors, we have not seen any sign 
of finite-volume correction in any of these low moments. Furthermore the results for momentum, 
{ x )u-cU an d helicity, (jt)A«-Arf> fractions obtained from a smaller, 1.8-fm, box at the same cutoff, 
agree with the larger box results. And at the lightest pion mass of 330 MeV these two quantities 
show interesting trending toward the experiments. Thus it would be interesting to study these 
moments at lighter mass. 



4. Iwasaki+DSDR ensembles 

We have been generating a new series of (2+l)-flavor dynamical DWF ensembles in a larger 
physical volume of about (4.6fm) 3 [9]. This was made possible by moving to a coarser lattice 
cutoff of cT x = 1.368(7) GeV. If we were to use such a coarse cutoff with Iwasaki gauge action 
we used for earlier studies, then we would have suffered from too large residual violation of the 
chiral symmetry. Thus we changed our gauge action by multiplying it with a ratio of Wilson- Dirac 
fermion determinants [13, 14, 15] that adequately suppress gauge dislocation that would result in 
unacceptably large residual violation of chiral symmetry while maintaining total gauge homotopy 
distribution comfortable for ergodicity of Monte Carlo gauge update algorithm. 

Two ensembles are being generated, one with degenerate up and down quark mass of 0.001 
lattice units and another at 0.0042. The strange quark mass is set at almost exactly physical point, 
at 0.0045 lattice units. With the observed residual mass of about 0.002 lattice units, the former 
ensemble corresponds to about 180 MeV pion mass, and the latter to about 250 MeV. We have 
accumulated about 1,100 hybrid Monte Carlo time units for the former and almost 1,900 for the 
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Figure 1: Preliminary nucleon effective mass with Gaussian smeared source of width 4 and 6 lattice units 
and point sink: the light, m n = 180 MeV, (left) and heavy, 250 MeV, (right). The lighter-ensemble results 
are from about 30 configurations each, while the heavier is from about 40. 



latter. Considering the thermalization time of 500 for the former and 600 for the latter, and the 
observable calculation interval of 8 time units, there are about 75 usable configuration for the 
former and 160 for the latter. At the time of the conference we had analyzed about 30 from the 
former and 20 from the latter. The statistics for the former stays the same while we have increased 
to about 40 for the latter now. 



5. Lattice systematic errors 

From our earlier studies we have identified two important sources of systematic errors [12]: 
1) time separation between nucleon source at time t = and sink at t^, and 2) spatial volume. 
In the above we have already explained how the latter is a serious problem for axialvector-current 
form factors but not for the other nucleon observables we are interested in. For the axialvector- 
current form factors the lattice spatial extent of five or more pion Compton wave lengths would be 
sufficient. Since our lightest pion mass now is about 1 80 MeV and the Compton wave length is 
about 1.1 fm, the spatial extent of about 4.6 fm may not be sufficiently large. But it is sufficient for 
the heavier pion mass of about 250 MeV. Thus we should at least be able to study the finite- volume 
correction in more detail. 

The source-sink time separation concerns excited-state contamination. No matter how we 
choose the nucleon operator, it would create some excited-state components at the source, t = 0. 
Though such contamination should decay more quickly in time than the ground state, it can still 
contaminate our calculation if we set the sink too close to the source in time. In our RBC 2-fiavor 
dynamical DWF study [4] we have demonstrated that such excited-state contamination can be 
serious in some of the nucleon observables we are interested in: results for the isovector momentum 
fraction, {x) u -d shifted beyond statistical error when we changed the source-sink separation from 
about 1.1 fm to 1.3 fm. Of course the optimum separation depends on the details of the nucleon 
operator and the nucleon mass spectrum on the given gauge ensemble. So prior to the three-point 
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Figure 2: Nucleon mass, m^, plotted against pion mass squared, m^, from RBC+UKQCD (2+l)-fiavor dy- 
namical DWF ensembles. The latest and lightest two points at m K = 250 and 180 MeV, though preliminary, 
seem trending downward toward the experiment. 

correlator calculations it is imperative to optimize the combination of nucleon operator and source- 
sink separation. 

We chose to use a standard nucleon operator, N = e a bc{u^Cy^db)u c , and optimize its Gaussian 
smearing [18, 19]. For both the light, m n = 180 MeV and 250 MeV ensembles we have compared 
the Gaussian width of 4 and 6 lattice units. In Fig. 1 we present the nucleon effective mass with 
point sink obtained from the current preliminary statistics about 30 for the light and 40 for the 
heavy ensembles. In both ensembles the wider, width-6, results settle on a plateau more quickly, 
and the narrower width-4 results merge with them. From these observations we conclude the wider 
width-6 Gaussian smearing is sufficient for our study. 

With the present preliminary statistics of 30 or 40, the plateaux decay beyond 8 or 10 lattice 
spacings or about 1.1 or 1.4 fm. But since we have accumulated 75 configurations for the light 
and 160 for the heavy ensembles, and since we are adding two more source positions in time at 
t = 16 and 48 in addition to the current and 32, we will soon quadruple our statistics: source-sink 
separation of about 1.4 fm should be available. 

6. Summary 

We are starting our nucleon structure calculations using the new series of RBC/UKQCD joint 
(2+l)-flavor dynamical DWF ensembles [9] with pion mass as low as 180 and 250 MeV and lattice 
spatial extent as large as 4.6 fm. As is shown in Fig. 2, preliminary estimate for nucleon mass 
is 0.721(13) lattice units or 0.98 GeV when pion mass is 180 MeV, and 0.763(10) or 1.05 GeV 
when pion mass is 250 MeV. They seem trending down toward the experiment more strongly than 
our previous studies at heavier pion mass. We have optimized our Gaussian smearing of nucleon 
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source, and will be reporting results for isovector form factors and some low moments of isovector 
structure functions in the near future. 
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